Mixed-state evolution in the presence of gain and loss 
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A model is proposed that describes the evolution of a mixed state of a quantum system for which 
gain and loss of energy or amplitude are present. Properties of the model are worked out in detail. 
In particular, invariant subspaces of the space of density matrices corresponding to the fixed points 
of the dynamics are identified, and the existence of a transition between the phase in which gain 
and loss are balanced and the phase in which this balance is lost is illustrated in terms of the time 
average of observables. The model is extended to include a noise term that results from a uniform 
random perturbation generated by white noise. Numerical studies of example systems show the 
emergence of equilibrium states that suppress the phase transition. 
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Over the past decade there have been considerable re- 
search interests, both theoretical and experimental, into 
the static and dynamic properties of classical and quan- 
tum systems for which gain and loss are present [TJ [5] ■ 
This is in part motivated by the realisation that when a 
system is placed in a configuration in which its energy or 
amplitude is transferred into its environment through one 
channel, but at the same time is amplified by the same 
amount through another channel, then the resulting dy- 
namics can exhibit features that are similar to those seen 
in Hamiltonian dynamical systems. The time evolution 
of such a system can be described by a Hamiltonian that 
is symmetric under a space-time reflection, that is, in- 
variant under the parity-time (PT) reversal. 

Interests in the theoretical study of PT symmetry was 
triggered by the surprising discovery that complex PT- 
symmetric quantum Hamiltonians can possess entirely 
real eigenvalues [3], One distinguishing feature of PT- 
symmetric quantum systems is the existence of phase 
transitions associated with the breakdown of the symme- 
try. That is, depending on the values of the matrix ele- 
ments of the Hamiltonian, its eigenstates may or may not 
be symmetric under the parity-time reversal. In the 'un- 
broken phase' where the eigenstates respect PT symme- 
try, the eigenvalues are real and there exists a similarity 
transformation that maps a local PT-symmetric Hamil- 
tonian into a typically nonlocal Hermitian Hamiltonian 
[H|5]; whereas in the 'broken phase' the eigenvalues con- 
stitute complex conjugate pairs. The transition between 
these phases is similar to second-order phase transitions 
in statistical mechanics, accompanied by singularities in 
the covariance matrix of the estimators for the param- 
eters in the Hamiltonian (which can be seen by using 
methods of information geometry [B]). 

Experimental realisations of these phenomena are mo- 
tivated in part by the observation that the presence of 
a loss, which traditionally has been viewed as an unde- 
sirable feature, can positively be manipulated so as to 
generate unexpected interesting effects (cf. [7]). In par- 



ticular, PT phase transitions have been predicted or ob- 
served in laboratory experiments for a range of systems; 
most notably in optical waveguides [8HTU]. but also in 
laser physics [TTMl3| . in electric circuits [21 [15], or in 
microwave cavity [16j . 

In the quantum context, theoretical investigations into 
the properties and predictions of the evolution equation 
where gain and loss are present have thus far been con- 
fined primarily to pure states. However, quantum sys- 
tems are more commonly described by mixed states, 
especially when they are subjected to manipulations. 
Therefore, in order to describe or predict the behaviour 
of such quantum systems it is necessary to understand 
how a given mixed state might evolve. The purpose of 
the present paper is to propose a model that describes 
the evolution of a generic density matrix in this context. 
Our model is given by the dynamical equation 

f t =-i[H,p}-([r,p} + -2tr( P r)p), (1) 

where H = W is the Hermitian part of the Hamilto- 
nian generating ambient unitary motion, r = is the 
skew-Hermitian part of the Hamiltonian governing gain 
and loss, and [r,p] + = Fp + pT denotes the symmet- 
ric product. We refer to ([I]) as the covariance equation, 
on account of the structure of the term involving r . In 
this paper we shall motivate equation ([!]), investigate its 
properties in detail, identify the structure of stationary 
states, and show that a PT phase transition manifests 
itself by means of the time average of observables. 

An experimental implementation of an evolution of a 
purely quantum system requires a careful balancing of 
loss and gain. It should be evident, however, that the 
realisation of such an evolution for a pure state is diffi- 
cult, to say the least, since energetic manipulations of a 
quantum system inadvertently perturb the system. To 
account for the possible impact of uncontrollable ambi- 
ent noise, we consider an extension of the model ([I]) that 
includes an additional term generated by uniform Gaus- 
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sian perturbations. We shall study dynamical aspects of 
the extended model, identify the existence of equilibrium 
states by means of numerical studies, and show how PT 
phase transitions can be suppressed by noise. 

Let us begin by motivating the introduction of the co- 
variance equation 0. In the case of pure states, the 
norm-preserving evolution equation generated by a com- 
plex Hamiltonian K = H — if is given by 



dt 



= -i(H-(H)M-(r-(r)M. 



(2) 



Together with an additional equation for the overall prob- 
ability N — —{r)N, equation is equivalent to the fa- 
miliar complexified Schrodinger equation. The benefit of 
([2]) is that it is defined on the projective Hilbert space. 
Besides the covariance equation 0, there are many alter- 
native evolution equations for mixed states that reduce 
to when restricted to a pure state satisfying p 2 = p. 
An example is given by the double-bracket equation 



dp 
dt 



-i[H, P ]-[[r )P ], P ]. 



(3) 



For r = H, the norm-preserving equations and 
have been considered by Gisin (T7] as candidate equations 
to describe dissipative quantum evolution. That and 
([3]) are identical for a pure state can be seen by observing 
that pTp = ti(pr)p when p 2 — p. Among the possible 
dynamical equations, the covariance equation is sin- 
gled out on account of the fact that its formal solution, 
given an initial state po, can be expressed in the form 



Pi 



e -i(H-ir)t poe i(H+ir)t 
tT(e- i ( H - i nt poe KH+ir)t 



(4) 



which provides a natural generalisation of its unitary 
counterpart when 7^ = 0. Similarly, the dynamical equa- 
tion satisfied by an observable (F) — tr(Fp t ) reads 



d(F) 
dt 



= i([H,F])-([r,F] + )+2(r}(F}, 



(5) 



which agrees with the complex extension of the Heisen- 
berg equation of motion obtained in Rcf. 18 20J for pure 
states. Note that the covariance-type structure in 
has also appeared in the contexts of approach to ther- 
mal equilibrium [5T], dissipative motion [12J[22], and con- 
strained quantum dynamics [2"3] . 

Key properties of the evolution equation can be 
summarised as follows: (i) It preserves the overall proba- 
bility so that tr(pt) = 1 for all t > 0. This can be checked 
by verifying the relation dttr(p) — 0. (ii) Unlike a uni- 
tary time evolution, it does not in general preserve the 
purity of the state. In particular, we have 



-trp 2 = -4(tr(r P 2 )-tr( P r)tr(p 2 )), 



dt 



(6) 



and in general the right side of is nonzero when p 2 7^ 
p. Thus, the purity of the initial state is not preserved 



by when p is not a fixed point of the dynamics, but 
an initially pure state will remain pure. When H = 0, 
we have the relation 



dt 



tv(pr) = -2vax p (r) < 0, 



(7) 



from which it follows that: (iii) The imaginary part 
r of the Hamiltonian drives every state towards the 
ground state of r. (iv) The evolution equation is 
'autonomous' in the sense that the dynamical trajectory 
in the space of density matrices is determined uniquely 
by the specification of the initial density matrix po, and 
is not dependent on the kind of probabilistic mixture the 
initial state might represent, (v) In the case of a unitary 
motion, the speed v = tv{d t ^) 2 of the evolution of the 
state is constant of motion and is given by the Wigner- 
Yanase skew information v = 2ti(H 2 p) — 2tr(Hy/pHyfp) 
[25] . which reduces to the Anandan-Aharonov relation 
V = 2AH 2 for pure states ^fp = p. When the dynamics 
is governed by a complex Hamiltonian H — i_T, the evo- 
lution speed is not a constant of motion and is given by 
the expression 

v = 2(tr(H 2 p)-tr(H^pH^p))-2itr([H,r]p) 

+2 (tr(r 2 p) + tr(r^pr^p) - 2(tr(/») 2 ) , (8) 

which reduces to v = 2AH 2 + 2Ar 2 - 2i([r, H}) for pure 
states. 

Next we identify the fixed points of the dynamics. We 
begin by the following elementary observation concerning 
eigenstates of a generic complex Hamiltonian K = H — 
ir. Suppose that \4>) is a normalised eigenstate of K such 
that 



K\<t>)=\\ct>). 



(9) 



Then, if A — E—ij, where E, 7 real, we have (<j>\H\(j)) = E 
and {(j)\r\4>) — 7. In particular, an eigenfunction \<p) of a 
complex Hamiltonian K has a real eigenvalue if and only 
if (<fi\r\(j>) = 0. With this in mind, we establish the fol- 
lowing results: (vi) Every eigenstate of the Hamiltonian 
K, irrespective of whether the eigenvalue is real, is a fixed 
point of the motion 0. Furthermore, they are the only 
stationary states that are pure, (vii) The mixed station- 
ary states of the dynamical equation consist of convex 
combinations of the eigenstates of the the Hamiltonian 
K associated with real eigenvalues. Furthermore, the to- 
tality of mixed stationary states lies on the subspace of 
density matrices for which tr(_T p) = 0. 

The statement (vi) can be established as follows: From 
Q, if we set p = \4>){4>\, where \4>) is an eigenstate of K 
with eigenvalue E — ry, then we have 



AKt p e iKH 



-2yt 



(10) 
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iKt Po e iKH ) 



and hence tr(e 
gether, it follows that p t 



Po 



e 2 ~ ft . Putting these to- 
To establish (vii) we set 



Po 



,Pk\<t>k){<l>k\ 



(11) 



k=l 



where {pk} are nonnegative numbers adding to unity, 
{10*;)} are eigenstates of K with real eigenvalues, and m 
is the number of real eigenvalues. Substituting (11) in 
Q, a short calculation shows that p t — pq. To show 
that the stationarity breaks down if the convex combina- 
tion contains eigenstates with complex eigenvalues, sup- 
pose that we add terms in |TT| ) associated with complex 
eigenvalues. Then from Q we obtain 



Pi 



-2 7j -t 



where -fj = (<f>j\r\<f>j)- Since at least one of the {jj} is 
nonzero by assumption, it follows that p t ^ po- Finally, 
since all terms in ( |TT| ) have the property that {<pk\r\<pk) = 
0, it follows that tr(Pp) = for all stationary mixed 
states. 

A corollary to (12) is that: (viii) If the initial state po 



admits an eigenfunction expansion such that one or more 
terms are associated with complex eigenvalues, then 



lim p t 

t—too 



WW 



(13) 



where |</>*) is the member of the eigenf unctions {|0j)} in 
the expansion for which the associated imaginary part of 
the eigenvalue jj takes maximum value. Hence the flow 
structure in the space of density matrices, when there are 
imaginary eigenvalues, is rather complex and intricate in 
higher dimensions, where for each eigenstate |0)(0| asso- 
ciated with an imaginary eigenvalue there is a continuum 
of states for which |0)(0| is an asymptotic attractor. If, 
however, there is no real eigenvalue at all, then this seg- 
mentation disappears and the eigenstate with the largest 
7 becomes the single attractor. 

To investigate the behaviour in the PT symmetric 
phase where all eigenvalues are real, let us put 



Po 



^Pjk\<i>j){<t>k 



(14) 



3,k 



Then the trace condition tr(po) = 1 implies that 

m 

j = l j^k 



(15) 



since the eigenfunctions are not orthogonal when f^O. 
Substituting (fl4j) in @ wc find 



Pt 



(16) 





(a) unbroken phase 



-0.5 °-5 

(b) broken phase 



FIG. 1: (colour online) Trajectories generated by the Hamiltonian 
K = a x — i-fUz- In the unbroken phase (|7| < 1) all orbits are 
closed, and the fixed-points are represented by a straight line. In 
the broken phase (I7I > 1) a sink and a source emerges and all 
states converge to the ground state in the limit t — > 00. 



(12) where uijk = Ej — E^. We thus deduce the following: (ix) 
When all eigenvalues of the Hamiltonian K are real, every 
orbit is periodic if the energy eigenvalues are commensu- 
rable; otherwise, the evolution is typically ergodic on a 
small toroidal subspace of the space of density matrices 
containing p Q . Therefore, in the PT-symmetric phase, 
dynamical features of the system are analogous to those 
of a unitary system, albeit differences in detail such as 
the non-constancy of the evolution speed. The similarity 
is due to the fact that the nonlinearity of ([!]) is merely 
to preserve tr(p); hence the evolution equation cannot 
generate nontrivial fixed points such as saddle points. 

In figure [I] we illustrate integral curves of ([!]) for a 
two-level example system with the Hamiltonian K = 
a x — i"/a z , where 7 is a real parameter. Note that the 
eigenvalues of K are given by A± = ±\/l — 7 2 - Hence 
the symmetry is broken if I7I > 1. Observe that in the 
unbroken phase the orbits are periodic, and hence the pu- 
rity of mixed states oscillates, while in the broken phase 
every mixed state is asymptotically purified to the ground 
state. To visualise the effect of the phase transition it will 
be useful to introduce an 'order parameter' m by the time 
average of a z : 



m = lim 



ti(a z p t )At. 



(17) 



Note that m is independent of the choice of the initial 
condition pq. In figure 2(a) we plot m as a function of 
7 , showing the characteristic behaviour of the order 
parameter in a second-order phase transition. 

Perhaps a surprising feature of the foregoing analysis 
is the observation that in the broken phase where eigen- 
values are complex, convex combinations of energy eigen- 
states are not stationary. On account of the fact that the 
dynamical equation is autonomous, this implies, in par- 
ticular, that if an initial state po is chosen randomly out 
the eigenstates \4>j)((f)j\, each one of which is stationary, 
and the system evolves under the presence of gain and 
loss such that PT symmetry is broken, then after a pas- 
sage of time the statistics obtained from p t is different 
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from those obtained from p$. This feature has no ana- 
logue in the standard unitary theory, and can be used 
to test the applicability of the model in laboratory 
experiments, provided that a coherent control of the sys- 
tem, in such a way that gain and loss can be balanced 
without perturbing the system, is possible. 

If, on the other hand, a coherent implementation of 
gain and loss is not feasible, either because of fundamen- 
tal quantum limits or current technological limits, then it 
is important to take into account additional effects aris- 
ing from random perturbations. For this purpose, we 
shall assume that the model remains valid, but in ad- 
dition the system is perturbed at random. Specifically, 
we assume that the state is perturbed by Gaussian white 
noise in every orthogonal direction in the space of pure 
states, with strength s/k, where k > 0. In Ref. [25] it 
was shown that such a perturbation, when averaged out, 
generates a flow in the space of density matrices given by 
k (1 — np); this, in turn, leads to the extended model: 

d/ ' - -i[H,p]-([r,p]+-2i X ( J >r)p) 



dt 



+k (1 — np) 



where k > and n is the Hilbert space dimensionality. 
We proceed to analyse properties of |l8|). To begin, it 



should be evident that: (i) The evolution equation (18) 



preserves the overall probability so that tr(p t ) = 1 for all 
t > 0. On account of the presence of noise, however, an 
initially pure state ceases to remain pure. In particular, 
we have: (ii) The evolution of the purity is governed by 
the equation 



dt 



tip 2 = -4(tr(Tp 2 ) - tr(pr)tr(p 2 )) 
+2k(1 - ntr(p 2 )). 



(19) 



Notice that when p 2 = p the first term in the right side 



of (191 vanishes, whereas the second term is negative. 



Hence an initially pure state necessarily evolves into a 
mixed state due to noise, (iii) When F = 0, the solution 



to the dynamical equation ( 18 ) takes the form 



Pt = -[l + (ne~ im p Q j m 



lie" 



(20) 



and has a single fixed point p^ — n~ 1. This can be seen 
from the facts that the off diagonal elements of n (1 — np) 
are negative, and that the diagonal elements are positive 
(reps, negative) if the diagonal element of p is less than 
(reps, larger than) The effect of the term k (1 — np) 
therefore is to generate a gradient flow towards the uni- 



formly mixed state p a 



L l. (iv) Like the model (h|, 



the evolution equation ( 18 ) is autonomous in the sense 



described above, (v) The dynamical equation satisfied 
by an observable (F) = tr(Fp t ) is given by 

d{F) 



dt 



i{[H,F])-{[r,F} + ) 
+n(tr(F) - n(F)). 



2(r)(F) 



(21) 




(a) re 



(b) re + 



FIG. 2: (colour online) Time and ensemble averages of the observ- 
able o~ z as functions of 7~ 1 . (a) When re = 0, irrespective of the 
initial state po, the time average of tr(o z pt) shows the existence of 
a phase transition at the critical point 7 C = 1. (b) When re ^ an 
equilibrium state p* emerges for each values of 7, re, to which every 
initial state po converges. Hence a time average can be replaced by 
an ensemble average. For re > the phase transition is eliminated, 
as indicated here for re = 0.01 (black), 0.1 (blue), 1 (magenta), and 
10 (green). 



To gain insights into fixed-point structures of the ex- 



(18) tended model (18) we have performed numerical stud- 



ies based on a two-level system with the Hamiltonian 
K = a x — v^u z - Our analysis shows that once the noise 
strength k is turned on, irrespective of its magnitude, the 
phase transition is suppressed. This can be seen by the 
consideration of the order parameter m. As indicated 
above, in the noise-free case there is a clear indication 
of a second-order phase transition at the critical point 
7 C = 1. Under a noisy environment, however, the sym- 
metry if broken for all values of 7 7^ 0; instead, for each 
value of 7, k an equilibrium state p* is established, to 



which all initial states converge. In figure 2(b) we plot 
m as a function of 7 _1 for a range of values for n > 0, 
showing the removal of the phase transition, although 
for sufficiently small k the signature of the transition is 
visible. We have performed further numerical and ana- 



lytical studies of the evolution equation (18), the results 



of which indicate that the emergence of a nontrivial equi- 
librium state is a generic feature of the model. 

In summary, we have introduced a model for describing 
the evolution of a density matrix for a system having gain 
and loss, and investigated its properties in detail. In 
particular, we have identified the associated stationary 
states, and shown the existence of phase transitions in 
the generic mixed-state context (cf. figure 2(a)[ ). We 
then extended the model to include a noise term, and 
showed evidences for the existence of equilibrium states 
that eliminate phase transitions (cf . figure |2(b)| . The 
applicability of our models are expected to be verifiable 
in laboratory experiments. 
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